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Abstract 

A kinetic theory is developed to describe radiating electrons whose motion is governed by the Lorentz- 
Dirac equation. This gives rise to a generalized Vlasov equation coupled to an equation for the evolution of 
the physical submanifold of phase space. The pathological solutions of the 1-particle theory may be removed 
by expanding the latter equation in powers of t := q'^ /&-Km. The radiation-induced change in entropy is 
explored, and its physical origin is discussed. As a simple demonstration of the theory, the radiative damping 
rate of longitudinal plasma waves is calculated. 

1 Introduction 

The question of how a charged particle interacts with its own radiation field remains unclear, despite intensive 
theoretical investigations over the last century. Until recently, this issue has been largely of theoretical interest as 
in current facilities the self-force is small in comparison to the Lorentz force due to the applied fields. However, 
the advent of new ultra-high intensity laser facilities requires that radiation reaction effects be taken seriously. 
For instance, the Extreme Light Infrastructure (ELI) [T] is expected to operate with intensities exceeding 10^^ 
W cm^^, and electron energies in the GeV range, at which level the radiation reaction force becomes comparable 
to and can even exceed the Lorentz force. 

The problem was first addressed by Lorentz P] and Abraham [5] , who introduced a third order differential 
equation to describe the trajectory of a radiating particle. This work was later generalized to the relativistic 
regime by Dirac [1], leading to what is now known as the Abraham-Lorentz-Dirac equation or, more simply, the 
Lorentz-Dirac equation, 

= --J^%(7'' + rA%c\ (1) 
m 

for the components {C°} of the worldline C : A i— )• = C"'(A) of an electron of charge q and mass m in an 
electromagnetic field Fab expressed in inertial coordinates {x"") on Minkowski spacetime M. with metric tensor 

7] = -dx° ® dx" + dx'^ (g) dx^ + dx^ (g) dx^ + dx^ ® dx^ . (2) 

T := /Qnm ~ 10^^'^s is the characteristic time of the electror0, A^b := 5^ + C°'Cb is the C-orthogonal 
projection, and an overdot indicates differentiation with respect to proper time A, 

JlabC'^C^ = -1, (3) 

where [r^ab] — diag( — 1, 1, 1, 1). We work in Heaviside-Lorentz units with c = 1, and raise and lower indices with 
rjab- Latin indices run from to 3, and the Einstein summation convention is used. 



♦Department of Physics, SUPA and University of Strathclyde, Glasgow, G4 ONG, UK 

^Department of Physics, Lancaster University, Lancaster, LAI 4YB, UK and Cockcroft Institute, Daresbury, WA4 4AD, UK 
^For definitencss we consider electrons, though the formalism could be applied to other charged particles. 



1 



Equation ([T|) has proved highly controversiaL The difficulties stem from the fact that it contains derivatives 
of the acceleration, so specifying the initial position and velocity is not sufficient to determine the solution. 
On the other hand, a generic specification of the initial acceleration leads to exponentially growing proper 
acceleration. 



C-Ca ^ e^l\ (4) 

even in the absence of applied forces. Such 'runaway solutions' are clearly at odds with our observations, and 
so must be eliminated. 

The solution proposed by Dirac was to use the elimination of the nonphysical solutions, rather than the 
initial acceleration, as part of the initial data, effectively converting the Lorentz-Dirac equation into a second 
order integro-differential equation, 

/•oo 

&{\)^ \ K''{X + aT)e-°'da, 
Jo 

K" := -—F'\C^ -T&CbC". (5) 
m 

Equation ([5]) requires the acceleration at a given time to depend on the force applied at all subsequent times, 
and so is acausal. However, the acausal effects are exponentially damped for times > r into the future, and so 
in practice are unobservable. 

Despite the many distinct derivations of the Lorentz-Dirac equation [fflllj , the problems of runaway solutions 
and acausality have led many researchers to propose alternative descriptions of radiation reaction |12H15) . The 
most widely adopted of these is the Landau-Lifshitz equation |16j . 



(7« = -^±F\C^ - -T\ddF\C'' - -A\F'',F^d\C^, (6) 
m m m 



obtained by perturbing ([T]) about the Lorentz force and neglecting terms higher than first order in r. It has 
been argued by Spohn [17| that the solutions approximated by (jG]) are precisely those corresponding to Dirac's 
asymptotic condition ([5]). Further, it is often claimed [TS] that effects neglected in the derivation of the LD 
equation - such as spin and quantum effects - are of order r^, so that ([5]) is as accurate as ([T|). Though 
plausible in the case of a single radiating electron, this claim remains unproven, and it is more reasonable to 
consider ^ an approximation to ([1]). There has been some interest recently in exploring the validity of this 
approximation [T91I20] . 

In practice, radiation reaction is unlikely to ever be observed in the context of a single radiating particle. 
However, modern laser facilities accelerate electron bunches with charge of the order of 10 pC, containing 10^ 
particles, and it follows that an appropriate description of radiation reaction is not the one-particle equation 
of motion, but rather a kinetic theory. A kinetic theory based on © has appeared in the literature pTVl23j . 
However, it is not clear a priori that the approximations leading to © are appropriate for describing a large 
number of interacting electrons. 

The Vlasov equation describes a large collection of particles as a continuum. Charged continua typically do 
not suffer from the pathologies associated with radiating point particles, and since the Maxwell- Vlasov system 
satisfies the energy and momentum conservation laws that led to ([1]) it may seem that it docs not need modifying 
to describe radiation reaction. However, the continuum comprises 'particles' with infinitesimal mass and charge, 
{q,m) — > with q/m fixed, in which limit ([T]) becomes the usual Lorentz force equation. To describe a collection 
of real electrons with finite q and jti, the Vlasov equation must be modified so that elements of the continuum 
follow trajectories of the Lorentz-Dirac equation ([IJ. Such a kinetic equation has been introduced in |24| : 
however, in that work the physical phase space was not identified, so the 1-particle distribution is required to 
be distributional (in the sense of Schwartz). As well as impeding the physical interpretation, this has led to a 
number of errors, for example in calculating the rate of change of entropy. 

In this article, wc introduce a kinetic theory for electrons obeying ([l]). with the 1-particle distribution a 
regular function on the physical phase space, and explore some of its consequences. 
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2 Kinetic theory with radiation reaction 



The usual approach to obtaining an equation for a large collection of particles from the equation of motion for 
a single particle is to consider a 1-particle distribution on a subspace of TM . This distribution is then taken to 
be constant along the lifts of particle orbits to TA4 , which may be obtained as integral curves of a vector field 
on TM . This prescription is not quite general enough to incorporate radiation reaction, since the Lorentz-Dirac 
equation is third order. The approach adopted in the present article is to consider a 1-particle distribution on 
a subspace of the bundle with total space 

TM(STM^ U TpM © TpM (7) 

and whose fibres are two copies of the tangent space, one representing 4-velocity and the other representing 
4-acceleration, with (i'^) and the induced coordinates on the first and second copies of the tangent space. 
We are interested in the physical subspace Q C TA4 ® TAi given by 

Q = {{x,x,x) G TMS)TM\ ipi =0, ip2 = 0,i° > 0} (8) 

where 

^l^^iVabX^'x' + 1), (9) 

^2 = Vabi'^i^ (10) 

Worldlines parametrized by proper time satisfy Lpi — 0, from which it follows that their 4-velocity and 4- 
acceleration should be orthogonal as encoded in ip2 ~ 0. It is straightforward to show that solutions to the 
Lorentz-Dirac equation are integral curves of the vector field L, 



^ - •^"Tnr + ^"tTIT + t*'*"^" + ^"'(*" + -F'bi')] 7F^' (11) 
on TM®TM. 

Since L is tangent to Q {Lipi ~ and Lip2 = '2riabX°'x''tpi + T^^ip2) it can be expressed as the push-forward 
of a vector field Lq on Q, L = i.*Lq, where i : Q ^ TA4 TA4 is the inclusion map. The Liouville vector field 
Lq is calculated explicitly below. 

The 1-particle distribution must satisfy the Vlasov equation, which states that the 1-particle distribution 
is preserved under the flow of the Liouville vector field. However, it is important to note that it is not the 
particle density / which must be preserved along the flow, but rather the particle distribution fuj, where lo is 
a non-vanishing top-dimensional form. Thus, the Vlasov equation may be cast as 

/:ls(.M = o (12) 

where Cx is the Lie derivative with respect to X, / is a 0-form on Q, and a; is a non-vanishing 10-form on Q. 

In principle, lj may be any measure on Q; however, in practice, it is convenient to use the Leray measure of 
Q C TM ® TM derived from the natural measure lo on TM (S TM , 

^ = d2;"123^ ^^0123 ^^^0123^ (^3) 

where da;*^^^"^ = A dx^ A dx^ A dx^, etc. The Leray measure w is defined by 

uj L*Cb, (14) 
where a) is any 10-form on TM ® TM satisfying 

a) = tJ A dipi A dip2. (15) 
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We choose 

Cj = {X^i° - i0il)-ldx"123 ^ ^^0123 ^ ^^23^ (^g) 



(17) 



which can readily be shown to satisfy ([T5|) . 

To calculate we use coordinates {x°',v^^,a'^) on Q such that 



* ■ 

L X 



= VT+V^, L*X^' = W^, (18) 

where Greek indices run from 1 to 3 and a ■ v = a^^v^, etc. 

Using dH]), ((ni), (HH), (HH) it foUows that the Leray measure a; = (*c;; on Q is 

cj = ^-dx0i23 A dai23 A dv^^\ (20) 

Similarly, on the fibre Qp = 7r~^(p) of the bundle {Q,TT,Ai) over p e A^. the induced measure ft is identified as 
^ dai23 A d„123 

^ - ^ 5 21 

because we adopt dx'^^^^ as the measure on 
It may be shown that 

_d_ _ _d_ 

'*dx-~dx-' ^ ' 

d d V,, d ay f ^„ v^Vu \ d 



* dv^ (9if yrr^ dx° vTT^ i+vydio' ^ ^ 

, ^ - ^ I A (24) 

and it follows t*Lg = i where 

-'^^ + + (^^"^^"^ + ^"("'^ + ^^^^*")) (2^) 

In (pSj) . and from now on, it is to be understood that i" and x° are shorthand for the values of and 
respectively. 

In the absence of radiation reaction, the Leray measure on the total space £ is preserved by the flow of the 
Liouville vector field induced from the Lorentz force equation. This well-known result is commonly described 
as "conservation of phase space volume". By contrast, the Leray measure a; on Q is not preserved by the flow 
of Lq: 

3 

Clq^ = (26) 

T 

This may be understood physically as a consequence of losses due to radiation. It follows from that ([T^ 
can be written 

LqJ + -/ = 0. (27) 

T 



^The same procedure leads to dx^^^'^ A dv^'^'^ + v'^ for the Leray measure on the bundle with total space £ := {{x,x) S 
TM\ tfi =0,x° > 0}, where is the pull-back of if* to £ with respect to the inclusion £ ^ TA4. It follows that the induced 
measure on a fibre of that bundle is dv^^^ / vl + v^ as expected. 
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The electrons couple to the electromagnetic field through the final term in ((25|) . and through the Maxwell 
equations 

dx°- dx^ dx'^ 

+ -^e-, (28) 

where the electron current is given by 

■J" =<lj fi"^ (29) 
and Jgxt is the current of any source other than the electrons. 

3 Physical solutions 

Equation ([27]) does not in general reduce to the usual Vlasov equation in the limit r — > 0. This is because, for 
a regular solution to (j27p . / will be nonzero over a range of a for given values of {x, v), and so will necessarily 
include runaway solutions as well as physical ones. To avoid this, look for solutions of the form 



f{x, V, a) = Vl + v^g{x,v)6^^\a- A{x,v)), (30) 

where the (as yet undetermined) functions A^^ describe the submanifold of phase space containing physical 
trajectories. The factor y/l + v^ in ([50)) ensures that g can be interpreted as a particle distribution function on 
the space of unit timelike 4- vectors, from which it follows that the current ()29p becomes 



Using ([50)) in ([?f| and integrating out the acceleration variables leads to the coupled equations 
dAf^ dAf^ 1 

+ - ^''^"-^ + + ^^"'^i")' (32) 

i"|^ + v/TT^^f-l^)=0, (33) 
dx"' ^v^'\y/T^^J ^ ' 



with A" v^Af^/\/l + v^. Equation ([5^ describes the evolution of the physical submanifold, as governed 
by ([T|). while ([55)) is a generalized Vlasov equation for the 1-particle distribution g. It is the reduced system 
([32)) -([33 )) . rather than ([27)) . that reduces to the usual Vlasov equation in the limit t ^ 0. 

So far we have placed no restrictions on the submanifold represented by A'^, other than that it satisfy ([32)) . 
In order that it represent the physical solutions to ([1]), it must be regular in the limit t ^ 0. A^ may then be 
expanded in powers of t: 

oo 

^" = E-"^(n)' (34) 

n=0 

with 

dA^ ri n 

J=0 J=0 
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Although T appears explicitly only in the acceleration equation (j32|) . the distribution g and electromagnetic field 
Fab acquire a r dependence via their couplings to in the Vlasov equation ((33)) and the Maxwell equations 
(|28p . respectively, and may be similarly expanded: 

cjo oo 

.9 = E^"5(„), F'^'' = ^r"^^(t). (37) 

n=0 71=0 

Truncating at n = yields the usual Vlasov equation without radiation reaction, while truncating at = 1 
yields the kinetic theory derived from the Landau-Lifshitz equation ([5]). This supports the large body of 
work pTH23] utilizing the latter approach. 



4 Entropy 

In the absence of radiation reaction, the entropy of the electrons is conserved. This is not generally true if 
radiation reaction is included, as the velocity-divergence of the physical acceleration acts as a source of entropy. 
Consider the entropy density s ~ —glng on phase space, which from the Vlasov equation p3p satisfies 

■a 9s r d ( sA^ \ r 9 / \ 



Then the entropy current s° = / x"' sci?" v j \/ \ + on spacetime satisfies 

girri^^Wv, (39) 



and the change in the total entropy 5* = / s'^dx of the electrons is 

9l^.{^^^]d'yd'^. (40) 



Although the entropy density s is defined with respect to a reference density, conservation of particle number 
ensures that changing this reference density merely shifts the total entropy by a constant, so the rate of change 
of entropy PO]) is well-defined. 

Restricting to the physical solutions of ([5^ . we can expand S — X^^o '''"•^(ti)' ^^'^ follows that 



J = ^ V -L 

The contribution from A^'g^ vanishes, so the leading order entropy change comes from n = 1, 

^ = -- / (ja{r + J:.,)+^^,TabS^'')d'x, (42) 
at m J \ / 

where 

Tab = FacFb'' - \TlabFcdF''' (43) 

is the strcss-energy-momentum tensor of the electromagnetic field and 

S''^ = m ( g.T^i''— £^ (44) 

is that of the electrons. Note that J" and S""^ in (05]) are calculated in the limit r — > 0. 

The rate of change of entropy ((42|) has previously been derived in [21] from a kinetic theory based on 
the Landau-Lifshitz equation, though the physical interpretation was there less apparent. The self- interaction 
term, JaJ", increases the entropy, while the interaction with the field, TabS""^, leads to an entropy decrease. 
The effect of the remaining term, JaJ^^t, depends on the nature of any external currents. In [53], only the Jq J° 
contribution was found, leading to the incorrect statement that the entropy of the electrons always increases. 
In fact, it is quite possible for the TabS""^ term to dominate, leading to radiation cooling of the electrons. 
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5 Plasma oscillations 



A simple demonstration of the theory comes from exploring the effect of radiation reaction on plasma waves. 
Linearize ([25]) . about the background 

.9 = A^=0, Fab = 0, (45) 

with J°xj = —quoSq and no ~ j gcPv. Assuming longitudinal perturbations of the form expi(fcz — ujt) yields 
the dispersion relation 



where A = 1 + zt(w\/1 + ~ kv^) represents the modification of the plasma waves due to radiation reaction. 

The dispersion relation (j46p possesses a more complicated mode structure than its counterpart in the 
Maxwell- Vlasov system. As well as small modifications to the usual solutions, there exist entirely new roots. 
However, the latter do not exist in the limit r — > 0, and so must be rejected as unphysical. 

Taking the cold equilibrium 'g{v) = noS^^\v) yields 

Lol^^uj^il + iTUj), (47) 



where u>p — \J (p-nojm, is the plasma frequency. Though rather simplistic, this nicely illustrates the effects of 
incorporating radiation reaction into the Vlasov equation. Two (physically equivalent) roots of (|47)) . 



C^«±C^p(l--(a.pT)2)--w2^, (48) 

indicate a damping rate of order r and an order frequency downshift. There also exists a third (purely 
imaginary) root, 

WW-, (49) 

T 

which represents extremely rapid growth, without oscillation. This mode, which does not exist in the limit t — > 0, 
corresponds to the runaway solutions of the Lorentz-Dirac equation, and must be discarded as nonphysical. 
Results equivalent to (|3S]) may be obtained by truncating the expansion in t at n = 2, while the root (|49p is 
excluded in this approach. Similar results may be found for other types of waves in plasmas. 



6 Conclusions 

In summary, with the advent of ultra-high intensity laser facilities, such as ELI, it is important to have a reliable 
kinetic theory of radiating particles, incorporating radiation reaction. We have developed such a theory based 
on the full Lorentz-Dirac equation, and found that it reduces to the usual Vlasov theory and to a kinetic theory 
based on the Landau-Lifshitz equation in appropriate limits. As simple demonstrations of the theory, we have 
explored the effects of radiation reaction on entropy and on longitudinal plasma waves. 
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